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ABSTRACT 

This  is  a  continuation  of  the  work  of  assembling  and  testing 
the  equipment  necessary  to  conduct  gamma-gamma  angular  correlation 
experiments 3  started  by  Fort  A.  Verser,  Jr.  at  the  US  Naval  Post- 
graduate School  (1). 

No  substantial  changes  were  made  in  the  equipment,,  but  rather 
a  rerun  of  the  Cobalt-60  experiment  was  made  to  confirm  the  results 
obtained  in  the  previous  work,   A  comparison  of  the  previous  work 
confirms  that  there  is  a  definite  distortion,  in  the  shape  of  the 
observed  correlation  function. 

The  statistical  analysis  of  the  data  and  the  solid  angle 
corrections  required  for  evaluation  of  the  results  have  been  pro- 
grammed for  calculation  by  use  of  the  1604  Digital  Computer.   This 
paper  outlines  these  programs  in  detail. 

The  author  wishes  to  express  his  appreciation  for  the  assist- 
ance given  him  by  Professor  Harry  E.  Handler  and  for  his  guidance 
throughout  this  investigation.   Appreciation  is  expressed  to 
Professor  Edmund  A.  Milne  for  his  advice  and  comments.   Captain  T. 
R.  Abernathys  USMC,  was  very  helpful  in  providing  valuable  assist- 
ance and  discussions  on  various  aspects  of  the  computer  programing. 
The  cooperation  and  assistance  of  Mr.  Mervyn  C.  Brillhart  in  main- 
taining and  testing  the  electronic  equipment  was  indispensable. 
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CHAPTER  I 
INTRODUCTION 

The  decay  of  Cobalt-60  to  an  excited  state  of  Nickel-60  and 
the  angular  correlation  of  the  resulting  gamma -gamma  cascade  have 
been  studied  extensively  experimentally  (2,6,8393  anc*  others)  and 
the  results  of  the  correlation  measurements  agree  very  well  with 
the  theoretical  correlation  for  the  decay  scheme  shown  in  Figure  1. 
Consequently,  the  decay  of  Cobalt=60  is  now  used  as  one  standard 
for  the  testing  of  correlation  equipment. 

Angular  correlation  equipment  was  assembled  at  the  US  Naval 
Postgraduate  School  by  Verser  in  1960  for  the  investigation  of  the 
properties  of  low- lying  excited  nuclear  states.   His  testing  of  the 
operating  characteristics  of  the  equipment  by  the  angular  corre- 
lation of  the  cascading  gamma  rays  resulting  from  the  decay  of 
Cobalt-60  yielded  an  observed  correlation  curve  whose  shape  was  not 
in  agreement  with  the  theoretical  curve.   This  work  was  initiated 
in  order  to  attempt  to  verify  Verser 's  results. 

In  order  to  facilitate  the  statistical  analysis  of  the  data, 
a  program  has  been  worked  out  for  the  use  of  the  1604  Digital 
Computer.   The  program  follows  the  procedures  outlined  by  M.E. 
Rose  (2)  and  is  contained  completely  in  Appendix  1. 

An  additional  program  has  been  written  to  compute  the  solid 
angle  corrections,  which  must  be  applied  to  the  theoretical  curve, 
by  the  method  of  M.E,  Rose  (2).   This  progam  is  outlined  in  Ap- 
pendix II. 

The  results  of  this  work  agreed  within  statistical  expectation 
with  those  obtained  by  Verser,  indicating  that  the  unexpected  shape 
of  the  curve  is  apparently  real  and  not  an  instrumental  effect. 
These  results  are  given  in  Chapter  VI. 
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CHAPTER  II 
THEORETICAL  CONSIDERATIONS 

The  determination  of  the  angular  momentum  quantum  numbers s 
or  spins,  of  low-lying  short-lived  excited  nuclear  states  by 
the  nuclear  spectroscopic  tool,  known  as  angular  correlation  of 
successively  emitted  radiations,  is  quite  common  practice  today. 
Only  within  recent  years ,  because  of  the  vast  advancement  of 
electronics,  has  practical  utilization  of  this  tool  become 
possible. 

Since  nuclei,  under  ordinary  conditions,  are  randomly  ori- 
ented in  space,  it  is  not  possible  to  observe  a  radiation  pattern. 
This  is  due  to  the  fact  that  the  probability  of  emission  of  a 
radiation  by  an  excited  nucleus,  depends  on  the  angle  between  the 
nuclear  spin  axis  and  the  direction  of  the  emission.   Therefore, 
in  order  to  observe  a  radiation  pattern,  it  is  first  necessary  to 
orient  the  nuclei. 

One  method  of  orienting  certain  nuclei  consists  of  placing  the 
sample  in  an  electric  field  gradient  or  a  strong  magnetic  field  at 
a  very  low  temperature.   Another  method,  the  one  used  in  the  present 
consideration,  consists  of  selecting  nuclei  whose  nuclear  states 
have  spins  lying  in  preferred  directions.   This  happens  to  be  the 
case  for  the  intermediate  state  if  a  nucleus  decays  by  successive 
emission  of  two  radiations. 

In  an  angular  correlation  experiment  the  first  radiation  is 
observed  in  a  fixed  direction.   This  establishes  an  axis  to  which 
the  direction  of  emission  of  the  second  radiation,  originating  from 
the  state  formed  by  the  first,  can  be  referred.   The  second  radiation 
has  preferred  angles  of  emission  with  respect  to  the  direction  of 
the  first. 

In  order  that  the  correlation  exhibit  maximum  anisotropy,  the 
angular  momentum  vector  of  the  intermediate  state  must  not  change 


direction  significantly  before  the  second  radiation  occurs.   The  mean 

—  ft 
life  of  the  intermediate  state  must  be  less  than  about  10   seconds , 

the  typical  precession  period  of  the  nuclear  spin  about  local  per- 
turbing fields.   The  angular  correlation  function  is  determined  by 
the  nature  of  the  radiations  and  by  the  spins  of  the  states  involved, 
and  thus  its  measurement  may  lead  to  an  unambiguous  choice  of  spins 
for  the  states. 

The  theoretical  expressions  for  the  correlation  functions  have 
been  worked  out  for  a  large  number  of  cases  of  interest  (7) .   For 
the  particular  case  of  a  gamma-gamma  cascade,  the  correlation 
function  W(9)  is 


n  n 

max  max 


W(9)      _,  A  P  (Cos  9)   /. ,    B  Cosn9  (1) 

n  =  0  n  =  0 

where  P   is  the  Legendre  Polynomial  of  even  order  n,  n    is  the 
n  max 

smallest  even  integer  of  the  set  of  three  numbers  consisting  of 
twice  the  multipole  order  of  each  of  the  two  gamma  rays,  and  9  is 
the  angle  between  the  directions  of  emission  of  the  two  gamma  rays. 

The  theoretical  correlation  function  for  the  cascade  in  Nickel-60, 
assuming  point  detectors  and  a  point  source,  is 


W(9)  =  1  +  0.1020  P2(Cos  9)  +  0.0091  P4(Cos  9)         (2) 


with  an  anisotropy 


,Ov       „r„~0 


R  =  "(180  )  -  W(90  )    =  0.1667  (3) 


W(90°) 


Many  workers,  after  correcting  these  theoretical  results  for  the 
effect  of  non-point  detectors  have  found  good  agreement  between  them 
and  their  experimental  correlation  results. 


CHAPTER  III 
DESCRIPTION  OF  EQUIPMENT 

The  equipment  needed  for  the  gamma -gamma  correlation 
measurements  must  furnish  means  for  detecting  the  radiation, 
pulse-height  amplification  with  selection  or  discrimination, 
coincidence  analysis,  and  adequate  recording  of  the  counts  in 
the  various  counting  channels.   A  mechanical  system  must  also  be 
provided  for  mounting  the  detectors  and  changing  their  relative 
positions  and  for  holding  and  positioning  the  radioactive  source. 

The  apparatus  used  for  this  work  was  that  assembled  by  Verser 
and  described  in  detail  in  his  report.   Therefore,  only  a  brief 
description  of  the  equipment  and  of  the  minor  alterations  per- 
formed on  it  will  be  presented. 

Figure  2  is  a  photograph  of  the  mechanical  setup.   The  origi- 
nal base  of  the  correlation  table  has  been  replaced,  the  new  base 
is  of  solid  brass  with  a  bearing  mounted  shaft  to  enable  accurate 
manual  detector  positioning. 

The  electronic  equipment  (Figure  3)  has  the  functions  of 
detection  of  the  radiation,  pulse-height  amplification  and  dis- 
crimination, coincidence,  analysis,  and  recording  the  number  of 
counts.   The  only  change  in  the  electronic  equipment  was  in  the 
pulse-shaping  networks  used  in  the  initial  calibration  of  the 
apparatus.   The  pulse  shapers  were,  rebuilt  to  give  a  better  simu- 
lation between  the  test  pulses  and  the  actual  pulses  encountered 
under  operating  conditions. 
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CHAPTER  IV 
EXPERIMENTAL  ANALYSIS 

The  angular  correlation  function  W(9)  (Equation  1)  must  be 
expressed  in  terms  of  experimentally-observable  quantities.   The 
equipment  records  the  total  number  of  coincidences  as  a  function  of 
the  angle  between  the  detectors,  and  the  total  number  of  radiations 
captured  in  each  channel.   C  ,  the  total  number  of  coincidences  in 
a  counting  period,  is  given  by 


C  =  C  +  C  +  C,  (4) 

t    g    a    b 


where 


C   =  the  number  of  genuine  coincidences  due  to  the  two 
cascading  radiations, 

C  =  the  number  of  accidental  coincidences  due  to  the 
finite  resolving  time  of  the  electronic  equipment, 
and 

C,  -  the  number  of  background  coincidences. 


The  half -life  of  the  source  is  long  compared  to  the.  time  required 
to  collect  the  necessary  data;  therefore,  equation  (4)  can  be  re- 
arranged and  rewritten  in  terms  of  counting  rates; 


C  -  C  -  C   -  Cu  (5) 

g    t    a    b 

C  can  also  be  expressed  in  terms  of  other  parameters; 
g 

^  -   <£liaiieil  a22P22  +  ei2ai2P12£2ia  21*21> VWo^V    <6> 


where 


N   =  the  absolute  source  strength 
o 


the  fraction  of  nuclei  decaying  by  the  emission  of 
gamma  ray  "j"  and  gamma  ray  "k"  in  a  cascade 
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i  =  the  fractional  solid  angle  subtended  by  the  crystal 
in  Channel  "i" 

?.  .  =  the  fraction  of  gamma  rays  "j"  which  leave  the  source 
in  the  right  direction  to  enter  crystal  "i"  that  are 
not  absorbed  before  reaching  the  crystal 
£..  =  the  total  efficiency  of  the  crystal  in  Channel  "i" 
for  gamma  ray  "j" 

a. .  =  the  fraction  of  the  number  of  gamma  rays  "j"  captured 

by  the  crystal  in  Channel  "i"  whose  pulse  heights  are 

accepted  by  the  energy  discrimination  of  Channel  " i" 

W(9  )  =  the  angular  correlation  function  at  an  angle  9 
m  m 

averaged  over  the  solid  angles  used. 


N.,  the  counting  rate  in  Channel  "i",  can  also  be  expressed  as  a 
function  of  certain  of  these  same  parameters. 


where 


N.  -  A.  0.N  )       £.  .a.  .   P  .  .  +  N.,  (7) 

1      11  O' — r1      lj  lj    lj      lb 


0.  =  the  number  of  gamma  rays  "i"  per  disintegration 
N.,  -  the  background  counting  rate  in  channel  "i" 


From  Equations  (5)  and  (6)  we  get  for  the  correlation  function^ 

W(9  )  =  <—J={ r- >  (8) 

m  'fisfl      fN 

o  J  m  !     o    j  m 


where  f  is  the  appropriate  function  of  the  £'s,  a's,  ?'s,  and  0's 
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By  combining  Equations  (7)  and  (8),  we  get  for  the  correlation 

f unc  t ion 

/-.   .   . 

(C  -C  -C,  )  N  I 

tab    o 

W(0  )  =  {_, r ...... , _V  (9) 

<NrNlb>  <VN2b> 


where  g  is  a  function  of  the  £'s,  a's,  9!s,  A1  s  0!s« 

In  this  form  W(9  )  is  relatively  insensitive  to  equipment  stability. 
m* 

Therefore,  it  is  only  necessary  to  calculate  the  modified  correlation 
function  given  by 

r   W(9  )      [  (C  -C  -C  ) 

_  m      \   t   a   d 

g  o      /(NrNlb)  <N2-H2b)  /  ^ 

The  correlation  runs  were  done  with  the  movable  detector  at 
seven  angles:  90°,  105°,  120°,  135°,  150°,  165°,  and  180°.   A 
minimum  of  six  two-hour  counting  periods  were  completed  at  each 
angle,  and  the  angular  order  in  which  they  were  taken  was  varied 
to  minimize  any  systematic  fluctuations.   Each  run  yielded  single 
channel  counting  rates  of  approximately  1200  counts  per  second  and 
a  total  coincidence  rate  of  approximately  1.1  counts  per  second. 

Before  the  experimental  results  can  be  compared  to  theory,  the 
data  must  be  checked  and  corrected,  if  necessary,  for  the  various 
deviations  from  the  ideal  arrangement,  i.e.,  centered  point  source, 
no  spurious  coincidences  due  to  the  coincidence  resolving  time,  no 
scattering,,  no  background,  and  perfect  stability  of  the  electronic 
equipment.   In  addition,  the  theoretical  correlation  function  must 
be  corrected  for  non-zero  solid  angle  detectors.   The  methods  for 
accounting  for  these  effects  are  discussed  in  the  following  para- 
graphs. 


II 


A.  SOURCE   SIZE   AND   POSITIONING 

As  a  test  of  the  positioning  of  the  source  in  the  center  of  the. 
correlation  table,  the  integral  counting  rate  of  a  portion  of  the 
Cobalt-60  spectrum  was  measured  in  each  channel  as  a  function  of  the 
angle  between  the  counters.   This  was  done  by  Verser  and  checked  again 
after  the  base  of  the  correlation  table  had  been  changed.   Within  sta- 
tistical expectations,  the  counting  rates  were  constant.   The  sotirce  has 
a  volume  of  about  0.001  cubic  inches  and  can  be  regarded  as  a  point  source 
at  the  distance  at  which  the  detectors  were  placed  from  it . 

B .  SCATTERING 

Compton  scattering  occurring  outside  the  detectors  can  give  rise  to 
unwanted  coincidences.  These  in  general  will  tend  to  smear  out  the  mea- 
sured correlation  function.  This  problem  is  minimized  by  use  of  scintil- 
lation crystals  as  detectors  and  accepting  only  the  full-energy  peaks  of 
the  desired  gamma  rays.  The  method  for  determining  the  pulse  height  se- 
lection limits  has  been  discussed  in  detail  by  Verser. 

C.  BACKGROUND  CORRECTION 

With  the  discriminators  at  the  same  settings  as  were  used  during  the 
correlation  runs,  background  runs  were  made  before  and  after  the  corre- 
lation runs.   The  background  coincidence  counting  rate  was  negligibly  small, 
and  the  background  rate  of  each  channel  was  such  as  to  yield  total  back- 
ground counts  of  the  order  of  0.4  percent  of  the  total  counts  in  each 
channel  during  the  correlation  runs.   These  background  rates  were  used  as 
corrections  to  the  gross  counting  rates  obtained  during  correlation  runs. 

D.  DETERMINATION  OF  GENUINE  COINCIDENCE  RATE 

Since  the  background  coincidence  counting  rate  proved  experimentally 
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to  be  negligible,  Equation  (5)  becomes 


C   -  C   -  C 
g    t    a 


For  the  total  counting  time,  T,  the  total  number  of  coincidences  is  C  T 
and  the  standard  deviation  in  the  total  coincidence  rate  is 


°~t   ~ 


(11) 


Therefore,  the  standard  deviation  of  the  genuine  coincidence  rate  is 

C       C 

2     2     2      fc       a 

a-        =cr      +  <r      =     ___  +  ___  (12) 

g      t      a       x        T 

The  accidental  coincidence  rate  is  given  by 


C  =  21  N  L  '  21   N 
a    r   1  2    r   o 


(13) 


where 


X     =   Coincidence  resolving  time  of  the  coincidence  analyzer. 


In  order  to  minimize  the  accidental  coincidence  rate,  the  source 
strength  must  be  kept  as  small  as  possible.   However,  if  the  source 
strength  is  decreased,  the  counting  time  must  be  increased  to  obtain  a 
specified  precision  in  the  counting  rate.   Thuss  a.  compromise  must  be 
reached  between  precision  and  counting  time. 

The  smallest  resolving  time  obtainable  with  the  present  equipment  is 
experimentally  determined  to  be 


T  =  (0.1820  +  0.0016)  microseconds; 


this  resulted  in  an  accidental  coincidence  rate  that  was  about  30%  of  the 


13 


total  coincidence  rate.   The  approximate  source  strength  of  the  water 
solution  of  CoCl^  used  in  this  experiment  was  0.05  millicuries  of  Cobalt  60 

E.   REFERENCE  COUNTING  RATE 


With  the  discriminator  in  each  channel  adjusted  to  accept  only  the 
full-energy  peaks  of  the  desired  gamma  rays,  the  integral  counting  rate 
of  the  Cobalt-60  spectrum  in  each  channel  was  determined Q   This  gave  a 
reference  counting  rate  for  each  channel.   During  the  course  of  the  corre- 
lation runs,  either  the  discriminator  settings  or  the  gains  of  the  linear 
amplifiers  were  adjusted  before  each  run  so  that,  the  integral  counting 
rate  agreed  within  +  5  cps  with  these  reference  rates  (Table  1) . 

TABLE  1 


REFERENCE  COUNTING  RATES 


N  =  (1085.26  +  0.57)  cps 
N2  =  (1117.32  +  0,59)  cps 


As  a  measure  of  the  gain  fluctuations  during  each  counting  runs  the 
integral  counting  rate  of  each  channel  was  measured  after  the  run  using 
the  same  discriminator  setting  that  was  employed  during  the  run.   The 
difference  between  this  counting  rate  and  the  reference  counting  for  that 
channel  was  determined.   If  the  variation  exceeded  +_  20  cps  in  either 
channel ,  the  run  was  discarded. 

This  procedure  resulted  in  about  307o  of  the  runs  being  discarded 
because  they  fell  outside  of  the  maximum  acceptable  fluctuation.   This 
indicates  that  the  electronic  system  is  not  as  stable  as  it  should  be.   A 
check  of  the  line  voltage  revealed  substantial  fluctuations  during  any  24= 
hour  period.   The  temperature  of  the  room  was  not  controlled  in  any  manner 
and  varied  considerably.   These  could  have  been  the  major  causes  of  the 
large  counting  rate  fluctuations  since  the  discriminators  are  highly  sensi- 
tive to  both  temperature  and  voltage  variations. 
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F.   SOLID  ANGLE  CORRECTIONS 


The  solid  angle  correction  which  must  be  applied  to  the  theoretical 
correlation  function  to  account  for  the  non-zero  solid  angle  of  the  de- 
tectors and  enable  a  comparison  to  be  made  with  the  measured  function  is 
described  by  Rose  (2).   The  correction  for  two  cylindrical  crystals  whose 
axes  intersect  at  the  source  involves  the  numerical  evaluation  of  integrals 
of  the  form 


n,  1 


rx  (0) 

P   (Cos  0)   1  -e  Sin  0d0 

n  J 


(14) 


where 


T.  =  the  full-energy  absorption  coefficient  of  the  detector  in 

Channel  "i"  for  the  gamma  ray; 
^ .  =  the  half-angle  subtended  by  the  front  face  of  the  crystal 


in  Channel  "i" ; 

X.(0)  =  t.Sec0   for  O^0<  Tan 
i       x 


-1 


_r 


h.  +  t. 
l     l 


X.(0)  =  r.Csc0-h.Sec0   for  Tan' 


h.  i  t; 
i      i 


<0<^r. 


where 


h.  =  the  distance  from  the  source  to  the  crystal  in  Channel  "i": 

t.  =  the  thickness  of  the  crystal  in  Channel  "i", 
l 

r.  =  the  radius  of  the  crystal  in  Channel  "i". 
l 


The  attenuation  factors,  Q  ,  are  functions  of  these  integrals 


—         - 

-         -) 

I      , 

I          o 

Q     = 

n,l 

n,2 

n 

I      , 

I          o 

0,1 

os  2 

_     s    ... 

(15) 


The  corrected  correlation  function  becomes 

N 
max 


W(0)  =  \   '   Q  A  P  (CosG) 
X      .    n  n  n 


n=o 
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CHAPTER  V 
EXPERIMENTAL  RESULTS 

The  modified  correlation  function  for  the  k    run  at  angle  9»  given 

in   ° 

by  Equation    (10)    is 

,  C.    "    c      "    CK 
W.1     (0.  )    ■  a  b 

km 


(NrNlb)    (N2-N2b); 


This  means  that  the  seven  counting  rates  C  ,  C  ,  C,  ,  N,,  N1L ,  N~3  and  N„, 

t   a   b    1    lb   Z       zb 

must  be  determined  for  each  run  at  each  angle  and  the  standard  deviation 
of  each  rate  calculated. 

The  background  counting  rates  N-.,  and  N„,  for  each  run  were  experi- 
mentally determined.   The  single-channel  counting  rates  corrected  for 
scaling  losses  yielded  N,  and  N„.   The  accidental  counting  rates  were  cal- 
culated for  each  run  from  Equation  (13).   The  total  coincidence  counting 
rate,  C  ,  was  experimentally  determined.   The  background  coincidence  rate. 
C,  ,  was  negligible  throughout  the  experiment. 

The  experimental  value  of  the  modified  correlation  function  W.f  (9-) 

k  N  nr 

and  its  standard  deviation  were  calculated  for  each  run  at  each  angle. 

From  these  the  weighted  mean  W"  (9. )  and  its  standard  deviation  were  com- 

m 

puted  for  each  angle. 

The  weighted  mean  of  the.  measurements  taken  at  each  angle  was  deter- 
mined by  using  the  inverse  square  of  the  standard  deviation  of  each  mea- 
surement as  a  weighting  factor: 


N 


w 


I  (eJ 


W?eJ  -  -4? — f (18) 


< 


where 

N»  =  number  of  runs  at  the  angle  9' , 
m  m 
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The  standard  deviation  of  the  weighted  mean    (W")  was  calculated  from 
the  individual  standard  deviations: 


/ 


(fly)') 


i 


a; 


i 


A  comparison  between  the  standard  deviations  of  the  means  calculated 
from  Equation  (19)  and  those  calculated  from 


(T2K) 


•^ 


wk)-W(o 


Nm(Nm-/ 


(20) 


revealed  no  significant  differences.   Therefore  the  standard  deviation 
obtained  from  Equation  (19)  were  used  in  subsequent  calculations  and  the 
method  discussed  in  paragraph  E  Chapter  IV  for  treating  gain  instabilities 
is  realistic. 

The  values  of  the  weighted  means  and  their  standard  deviations  are 
tabulated  in  Table  2. 

TABLE  2 

The  Weighted  Means  W' (9. ) 

m 


Q. (degrees) 

Weighted 

Standard 

Mean 

Deviation 

90 

5046 

45 

105 

4925 

33 

120 

5109 

21 

135 

5211 

36 

150 

5556 

35 

165 

5667 

24 

180 

5907 

43 
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A  curve  of  the  form  W1 (9)  =  /    A1   P  (CosQ)  was  fitted  to  the  expert 


mental  W1 (9  )  by  the  method  of  least  squares  as  outlined  by  Rose  (2) 


m 


This  analysis  was  programmed  to  be  carried  out  by  the  1604  computer.   (See 
Appendix  II  for  details) . 

The  values  of  A'  are  given  by  the  matrix  equation 


T 

A  ojA 


1. 


II  atJ  •  IMI- 11/41 


(21) 


where 


=  the  7x3  matrix  consisting  of  the  Legendre  Polynomialsj 
P  ,  P9,  and  P, ,  evaluated  at  the  seven  angles  (Figure  5A) 

=  the  transposed  A  matrix 


j  u)||   =  a  diagonal  matrix  with  elements  consisting  of  the  inverse 

squares  of  the  standard  deviations  of  the  experimental 

W' (9,)  (Figure  5B) 
m 

fji\\     =   a  one  column  matrix  whose  elements  are  the  experimental 

values  of  W'(9.)  (Figure  5C) 

"  m 

A1!!   =  the  one  column  matrix  whose  elements  are  respectively 


-1 


the  least-square  coefficients  A* ,  Al  and  A!  (Figure  5D) 

A^tuAll   "*"=  the  inverse  of  the  matrix  obtained  by  multiplying  the 

three  designated  matrices  together.   The  diagonal  terms 

of  this  matrix  are  the  squares  of  the  standard  deviations 

of  the  least-square  coefficients.  A" . 

n 

The  experimental  values  for  the  coefficients  A8  and  their  standard 

n 

deviations,  normalized  so  that  A1  =  1,  are  listed  in  Table  3.   The  corre- 
5  o 

sponding  values  from  Verser's  work  are  included  for  comparison. 
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TABLE  3 
Least  Square  Coefficients 


Ver 

ser 

This  Work 

A' 
o 

1.0000  +  0.0034 

1.0000  +  0.0026 

A2 

0.0961  +  0.0055 

0.0971  +  0.0057 

Ai 

0.0339  +  0.0071 

0.0231  +  0.0066 
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MATRIX  CONFIGURATIONS 


P  (Cos90  ) 
o 

P  (CosI05°) 
o 

P  (Cosl20°) 
o 

P  (Cosl35°) 
o 

P  (Cosl50°) 

P  (Cosl65°) 
o 

P  (Cosl80°) 
o 


P2(Cos90  ) 

P  (Cosl05°) 

P  (Cosl20°) 

P2(Cosl35°) 

P2(Cosl50°) 

P2(Cosl65°) 

P2(Cosl80°) 


P4(Cos90  ) 

P  (Cosl05°) 

P  (Cosl20°) 

P4(Cosl35°) 

P4(Cosl50°) 

P4(Cosl65°) 

P4(Cosl80°) 


(A) 


(T 

'(90°) 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

cT"2(i05o 

)  o 

0 

0 

i 

0 

0 

(T 

2(120( 

3)  o 

0 

0 

0 

0 

0 

a- 

2(135°) 

0 

0 

0 

0 

0 

■ 

0 

C' 

2(150°) 

11 

0 

0 

0 

0 

0 

0"2(165°) 

0 

0 

0 

0 

0 

0 

tf"2(180O) 

(B) 


W' (90°) 
W' (105°) 
Wv (120°) 
W'(135°) 
W!(150°) 
W' (165°) 
W'(180°) 

(C) 


"4 
(D) 


FIGURE  5 
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A.   Experimental  Anisotropics 

The  anisotropy  is  defined: 


,ov    —,^o 


R  =  W(180  )  -  W(90  )  (22) 

W(90°) 

The  anisotropy  for  the  least-square  curve  was  calculated  from  the  experi- 
mental correlation  function.   Since  the  anisotropy  can  be  determined  with 
much  greater  precision  than  can  the  least-square  coefficients  themselves, 
it  provides  a  more  sensitive  means  for  comparing  the  experimental  and  then 
retical  functions.   The  experimental  anisotropy  is 

R  =  0.167  +  0.013 

The  standard  deviation  was  calculated  by  the  formula  presented  by  Kleroa 
and  McGowan  (6) .  The  computer  program  for  the  calculation  of  R  and  its 
standard  deviation  is  outlined  in  Appendix  I. 
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CHAPTER  VI 
RESULTS  AND  CONCLUSIONS 

A.   COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY 

The  experimental  and  corrected  theoretical  functions  for  this  experi- 
ment are: 

Experimental : 

W'(Q)  =  (1.0000  +  0.0026)  +  (0.0971  +  0.0057)  P? 
+  (0.0231  +  0.0066)  P4 


Theoretical : 

W(9)  -  1.0000  +  0„0996  P2  +  0.00838  ?^ 

The  experimental  anisotropy  was  calculated  to  be  R  =  0„167  +  0„013  and 
is  compared  with  Verser's  value  of  R  =  0.171  +  0.027  and  with  the  theo- 
retical anisotropy  corrected  for  the  solid  angle,  R  =  0„161. 

Figure  6  shows  the  experimental  results  and  their  standard  deviations  for 
this  work  and  that  of  Verser.   Also  shown  are  the  least-square  correlation 
curves  for  both  sets  of  data  and  the  corrected  theoretical  correlation 
curve.   All  the  curves  have  been  normalized  to  unity  at  0  =  90  . 

B.   DISCUSSION  AND  CONCLUSIONS 

From  Figure  6,  it  can  be  seen  that  the  curve  obtained  in  this  work 
agrees  with  that  obtained  by  Verser.   Both  of  these  curves  show  an  ap- 
parently real  dip  in  the  region  from  105   to  135  .   It  is  only  a  remote 
possibility  that  this  dip  can  be  attributed  to  statistical  fluctuations, 
and  it  is  unlikely  to  be  due  to  instrumental  effects. 
Two  conjectures  of  possible  mechanisms  are: 

(a)  The  source  may  have  crystallized  resulting  in  preferential 
alignment  of  the  nuclei  within  the  crystal  due  to  internal 
electric  fields. 
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(b)  The  ferromagnetic  properties  of  Cobalt  could  have  led  to 
a  preferential  alignment  of  the  nuclei  in  the  earth's  mag- 
netic field. 
The  following  series  of  experiments  are  suggested  for  future  work. 

(a)  The  results  for  measurements  conducted  for  angles  between 
180  and  270   should  be  compared  with  those  already  ob- 
tained. 

(b)  In  order  to  check  on  the  possible  effects  of  orientation 
due  to  crystalline  electric  fields,  the  experiment  should 
be  repeated  with  the  present  source  at  various  orientations 
relative  to  the  fixed  counter. 

(c)  In  order  to  check  on  the  possibility  of  orientation  by  mag- 
netic effects  the  source  should  be  placed  in  a  magnetic  field 
of  known  direction  and  the  measurements  repeated. 

(d)  A  correlation  curve  using  a  powdered  Cobalt-60  source  in 
place  of  the  present  source  should  be  obtained  for  comparison. 

These  experiments  would  determine  conclusively  if  any  effects  causing  pre 
ferential  nuclear  orientation  are  present. 

The  following  instrumentation  improvements  are  suggested  in  order  to 
improve  the  stability  of  the  equipment. 

(a)  A  very  stable  power  supply  should  be  purchased  to  aid  in 
preventing  the  line  voltage  fluctuations  from  interfering 
with  the  electronic  system. 

(b)  The  effect  of  temperature  variations  on  the  electronic  system 
can  be  minimized  by  placing  the  equipment  in  a  temperature 
controlled  room  for  conducting  future  wor 

(c)  A  slow  coincidence  unit  whose  resolving  time  is  independent 
of  the  counting  rate  should  be  obtained. 

The  precision  of  the  results  should  be  improved  by  the  addition  oi 
fast  coincidence  channel  working  in  conjunction  with  the  slow  coincidence 
unit  to  reduce  the  number  of  accidental  coincidences. 
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Since  the  experimental  curves  using  Cobalt-60  have  not  agreed  with 
the  theoretical  curve,  the  experiments  on  Cobalt-60  should  be  continued 
until  the  source  of  the  distrotion  is  ascertained.   Perhaps  the  effect 
can  be  related  to  the  physical  properties  of  the  CoCl9  source. 
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APPENDIX  I 
Least-Square  Fit  of  Experimental  Results 

The  comparison  of  the  experimental  correlation  function  to  the 

functional  form  of  the  theoretical  correlation  function  given  by  Equation 

(1)  is  accomplished  by  the  method  of  least-squares,  as  outlined  by  Rose 

(2) .   The  analysis  requires  the  solution  of  Equation  (20)  for  the  values 

of  A'  . 
n 

The  computer  solution  of  the  matrix  equation,  like  the  solid  angle 
calculation,  is  carried  out  in  floating  point  format.   In  addition  to  the 
solution  for  the  values  of  A' ,  the  program  also  calculates  the  anisotropy 
by  Equation  (21),  and  its  standard  deviation  by  the  formula  presented  by 
Klema  and  McGowan  (6).   A  flow  sheet  outline  of  the  program  is  given  in 
Figure  (7).   In  order  to  provide  a  certain  amount  of  flexibility  to  the 
program,  the  angles  used  have  been  left  as  a  parameter  to  be  supplied. 
The  program  as  written  is,  however,  limited  to  seven  different  angular 
positions o   The  parameters  needed  for  the  program  may  be  placed  anywhere 
in  the  computer  and  their  addresses  supplied  to  the  program  by  placing 
them  in  the  B  registers  in  the  following  order: 

Bl  =  weighted  means  of  the  correlation  function  in  ascending  angular 
order. 

B2  =  the  standard  deveations  of  the  weighted  means  in  the  same  order 
as  above. 

B3  =  the.  values  of  the  angles  used  in  ascending  order. 

The  program  is  available  for  future  use  on  both  bioctal  tape  in 
machine  language,  and  IBM  cards  in  assembly  language.   The  bioctal  tape 
contains  the  following: 

1.  The  main  program  for  calculation  of  A' 

n 

2.  The  program  of  normalization  of  the  A' 

3.  The  program  for  calculation  of  the  resulting  correlation  function 

4.  The  normalization  of  the  correlation  function 

5.  The  calculation  of  the  anisotropy 
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6„  The  calculation  of  the  standard  deviation  for  the  anisotropy 
7.    The  subroutines 

a.  Floating  Point  Sine-Cosine 

b.  Floating  Point  Square  Root 

c.  Floating  Point  Decimal  Output. 

The  results  of  the  computations  are  dumped  on  magnetic  tape  unit  (4) 

in  the  following  sequence: 

Lines  1  thru  3   The  matrix  inverse 

Lines  4  thru  6    Check  answer;  included  in  the  program  is  a  check 

of  the  inverse  which  consists  of  multiplying  the 

inverse  by  the  original  matrix 

Line  7  The  values  of  A' ,  A'  and  A' 

o    I  4 

Line  8  The  standard  deviations  of  the  A" 

n 

Line  9  thru  13   The  resulting  correlation  functions 

Line  14  and   15  The  normalized  correlation  functions 

Line  16  The  normalized  A1 

n 

Line  17  The  normalized  standard  deviations  of  the  A" 

n 

Line  18  The  anisotropy 

Line  19  The  square  of  the  standard  deviation  of  the  anisotropy 
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STATISTICAL  ANALYSIS  PROGRAM 


ORG 

12000 

REM 

BURTON  J  CONWAY  APRIL  I96I 

REM 

LEAST  SQUARE  FIT  OF  ANGULAR 

REM 

CORRELATION  DATA. 

REM 

MUST  BE  SUPPLIED  WITH 

REM 

PARAMETER  ADDRESSES  AS  FOLLOWS. 

REM 

Bl  CORRELATION  FUNCTION. 

REM 

B2  SIGMA  OF  CORRELATION  FUNGI*. 

REM 

B3  ANGLES  USED  IN  DEGREES. 

ZEUS 

EQU 

13000 

RHO 

EQU 

13100 

TRANS 

EQU 

13200 

ATW 

EQU 

13230 

ATWA 

EQU 

13260 

INVER 

EQU 

13500 

CHECK 

EQU 

13520 

IAT 

EQU 

135^0 

IATW 

EQU 

13570 

ANS 

EQU 

13620 

FINIS 

EQU 

13630 

SINFL 

EQU 

60000 

MATRIX 

EQU 

10000 
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12000 


12001 


12002 


12003 


12004 


12005 


12006 


1200? 


12010 


12011 


12012 


aniso 

EQU 

13624 

56  1  12045 

SIU 

1 

FORMMU 

50  i  ooooo 

ENI 

1 

0 

57  2  12052 

SIL 

2 

SIGMA 

56  2  12033 

SIU 

2 

SIGMA+1 

50  2  00000 

ENI 

2 

0 

50  0  00000 

ENI 

0 
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56  3  12004 
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3 
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50  3  ooooo 
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3 

0 

12  1  OOOOO 

START 

IDA 

1 

0 

32  0  12512 
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0 
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0 
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0 

0 

50  0  ooooo 

ENI 

0 

0 
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STA 

1 

TEMJ 

50  0  ooooo 

ENI 

c 

0 
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ISK 

1 

6 

75  0  12004 

SU 

0 

START 

12  1  12463 

NEXT 

IDA 

1 

TEMP 

STARTS  FORMATION 

32  1  12463 

FMU 

1 

TEMP 

OF  THE  LEG.  POL. 

30 


12013 


12014 


12015 


12016 


12017 


12020 


12021 


12022 


12023 


12024 


12025 


12026 


20  0  12461 

SEA 

0 

DON 

32  0  12320 

FMU 

0 

FLT2 

31  0  12316 

FSB 

0 

FLT 

33  0  12317 

FDV 

0 

FLT1 

20  1  12326 

STA 

1 

W0RK2 

12  0  12461 

LDA 

0 

don 

32  0  12322 

FMU 

0 

:-Li- 
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STA 

0 
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00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 

20  7  05015 
20  7  37137 


OCT 


OCT 


OCT 


OCT 


OCT 


CX  !T 


OCT 


OCT 


0 


0 


DEC    .4938271D+17 


OCT 


DEC    .4 52693 5D+17 
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1237^      00  0  00000  OCT    0 

00  0  00000 

12375  00  0  00000  OCT    0 
00  0  00000 

12376  00  0  00000         OCT    0 
00  0  00000 

12377  00  0  00000         OCT   0 
00  0  00000 

12U00     00  0  00000  OCT    0 

00  0  00000 
12^01      00  0  00000  OCT    0 

00  0  00000 
12U02     00  0  00000  OCT    0 

00  0  00000 
12U03     20  7  05176         DEC    .i+725897D+17 

27  1  0316U 
12U01+     00  0  00000  OCT    0 

00  0  00000 
121+05     00  0  00000  OCT    0 

00  0  00000 
12^06      00  0  00000  OCT    0 

00  0  00000 
12U07     00  0  00000  OCT    0 

00  0  00000 


51 


121+10  00  0  00000  OCT         0 

00  0  00000 
12U11      00  0  00000  OCT    0 

00  0  00000 

12412  00  0  00000  OCT    0 
00  0  00000 

12413  20  7  04477         DEC    .4l64931D+17 
37  2  35173 

12414  00  0  00000  OCT    0 
00  0  00000 

12415  00  0  00000  OCT    0 
00  0  00000 

12416  00  0  00000  OCT    0 
00  0  00000 

12417  00  0  00000  OCT    0 
00  0  00000 

12420  00  0  00000  OCT    0 
00  0  00000 

12421  00  0  00000  OCT    0 
00  0  00000 

12422  00  0  00000  OCT    0 
00  0  00000 
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12423  20  7  150*H 
l6  1  31752 

12424  00  0  00000 
00  0  00000 

12425  00  0  00000 
00  0  00000 

12426  00  0  00000 
00  0  00000 

12427  00  0  00000 
00  0  00000 

12430  00  0  00000 
00  0  00000 

12431  00  0  00000 
00  0  00000 

12432  00  0  00000 
00  0  00000 

12433  20  7  05367 
05  2  01171 

12434  20  0  14000   ABLE 
00  0  00000 

12435  57  7  73777 
77  7  77777 

12436  17  7  66000 
00  0  00000 


DEC    .9124087D+17 


OCT 


OCT 


OCT 


OCT 


OCT 


OCT 


OCT 


DEC 


DEC 


DEC 


0 


0 


0 


0 


DEC    .4938271D+17 


1.0 


-.5 


.375 


53 


12437 


12440 


12441 


12442 


12443 


12444 


12445 


L2446 


L2447 


12450 


L2451 


20  0  14000 
00  0  00000 
6o  o  11466 
73  1  10267 
17  7  54455 
71  7  14727 
20  0  14000 
00  0  00000 
60  0  23777 
77  7  77777 
60  0  13277 
76  5  40722 
20  0  14000 
00  0  00000 

17  7  64ooo 
10  7  23316 
60  0  11400 

05  1  74265 

20  0  14000 
00  0  00000 
20  0  05000 
00  0  00000 


DEC 


1.0 


DEC  -.  39954  51* 


DEC  .1434268 


DEC 


DEC 


J'C 


DEC 


DEC 


DEC 


DEC 


1.0 


,125 


DEC  -.289065 


1.0 


.250017 


-.40624 


1. 


.625 
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12452  17  7  26004 
31  o  11317 

12453  20  0  14000 
00  0  00000 

12454  20  0  07144 
26  4  1625^ 

12455  20  0  05364 
13  3  47506 

12456  20  0  14000 
00  0  00000 

12457  20  0  14000 
00  0  00000 

12460  20  0  14000 
00  0  00000 

12461  00  0  00000 
00  0  00000 

12463      00  0  00000   TEMP 

00  0  00000 
12472     00  0  00000 

00  0  00000 
12501     00  0  00000 

00  0  00000 
12504      00  0  00000   TTAN 

00  0  00000 


DEC    .023471 


DON 


NORM 


TTA 


DEC 


DEC 


DEC 


DEC 


DEC 


BSS 


BSS 


BSS 


BSS 


BSS 


1.0 


.8995 


DEC    .684614 


1.0 


1.0 


1.0 
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12507      00  0  00000   AN 
00  0  00000 

12512  1?  7  2U357   CONV 
50  6  50U51 

12513  00  0  00000   EBSIL 
00  0  00000 

1251^      00  0  00000   ROB 

00  0  00000 
12515      00  0  00000   AXE 

00  0  00000 

12520  00  0  00000   C0V02 
00  0  00000 

12521  00  0  00000   C0V2U 
00  0  00000 

12522  00  0  00000   A2U 
00  0  00000 

12523  00  0  00000  P2h2 
00  0  00000 

12524  00  0  00000   GET 
00  0  00000 

12525  00  0  00000   WIT 
00  0  00000 


BSS 


DEC    .01745329252 


BSS 


BSS 


BSS 


BSS 


BSS 


BSS 


BSS 


56 


12527 


12550 


L2531 


12532 


12533 


12531* 


12535 


12536 


12537 


125^0 


125^1 


75  o  ooooo 
56  1  12552 
22  0  12550 
22  3  12527 

50  o  ooooo 

U2  0  12556 
05  0  00001 
03  0  000**  5 

22  3  1253^ 
11  0  00002 

11  0  03777 
03  0  00001 

05  0  000UU 

20  0  I2561 

kk  0  12557 

01  0  0000U 
";.:  :  125U0 
01  0  00001 
20  0  12562 
i*+  0  12555 
20  0  12563 

12  0  1255*+ 


SQRFL 


REM 

FLOATING  P( 

)INT  SQUARE  ROOT 

SU 

0 

0 

ALARM  EXIT 

SIU 

1 

SQFEXIT 

STORE  Bl 

AJP 

0 

SQFEXIT-2 

TEST  0  ARGUMENT 

AJP 

3 

SQRFL 

TEST  -  ARGUMENT 

ENI 

0 

0 

SCM 

0 

SQMASKL 

SET  CORRECT 

ALS 

0 

1 

SIGN  FOR 

LRS 

0 

37D 

EXPONENT 

AJP 

3 

/+1 

COMPUTE 

INA 

0 

2 

EXPONENT 

INA 

0 

3777 

FOR  ROOT 

LRS 

0 

1 

ALS 

0 

36D 

SHIFT  EXPONENT 

STA 

0 

SQERAS 

AND  SIGN 

LDL 

0 

SQMASK2 

ARS 

0 

k 

QJP 

3 

/+1 

ARS 

0 

1 

STA 

0 

SQERAS+1 

ADD 

0  SQFC0N3 

STA 

0  SQERAS+2 

LDA 

0 

SQFC0N2 

57 


12542 


12543 


12544 


125^5 


12546 


12547 


12550 


12551 


12552 


12553 


12554 


12555 


27  o  12563 

DVF 

0 

SQERAS+2 

14  o  12553 

ADD 

0 

SQFC0N1 

50  1  00001 

ENI 

1 

1 

04  0  00000 

ENQ 

0 

0 

20  0  12563 

STA 

0 

SQERAS+2 

12  0  12562 

LDA 

0 

SQERAS+1 

27  0  12563 

DVF 

0 

SQERAS+2 

14  0  12563 

ADD 

0 

SQERAS+2 

01  0  00001 

ARS 

0 

1 

55  1  12544 

UP 

1 

/~2 

01  0  00011 

ARS 

0 

11 

42  0  12561 

SCM 

0 

SQJERAS 

32  0  12560 

FMU 

0 

SQHASK3 

52  1  12527 

LIU 

1 

SQRFL 

51  1  00001 

INI 

1 

1 

57  1  12552 

SIL 

1 

/+1 

50  1  00000 

SQFEXIT 

ENI 

1 

0 

75  0  00000 

SLJ 

0 

0 

21  3  66405 

SQFC0N1 

DEC 

218518306d=8b45 

01  4  47436 

76  3  72106 

SQFC0N2 

DEC 

-30228991714D-1 

43  6  47567 

03  0  56173 

SQFC0N3 

DEC 

154515776D-8B43 

52  7  71341 
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12557 


12560 


12561 


20  0  00000 
00  0  00000 
37  7  77777 
77  7  7^000 
20  0  14000 
00  0  00000 
00  0  00000 
00  0  00000 


SQMASK1  OCT 


2000000000000000 


SQMASK2  OCT   377777777777^000 


SQ4ASK3  OCT    2001^00000000000 


SQERAS   BSS 


END 
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APPENDIX  II 
Solid  Angle  Correction 

The  solid  angle  correction,  which  must  be  applied  to  the  theoretical 
correlation  function  to  enable  a  comparison  to  be  made  with  the.  experi- 
mental correlation  function,  is  described  by  Rose  (1) .   The  correction 
involves  the  numerical  evaluation  of  the  integrals  given  by  Equation  (14) 

Ini  =  (       Pn(Cos0)  (l-e"TiXi(0))  Sin0  d0 
•J  o 
and  the  subsequent  solutions  of  Equations  (15)  and  (16) . 

The  computer  solution  for  this  correction  was  carried  out  in  floating 
point  format.   The  numerical  integration  was  computed  using  an  interval 
of  0.001  radians.   A  flow  sheet  outline  of  the  program  is  given  in  Figure 
(8).   The  program  must  be  supplied  with  the  parameters  T »  h,  r,  and  t 
same  in  both  channels  (see  Table  4)  thus  enabling  solid  angle  corrections 
to  be  made  for  any  experimental  set-up  to  which  these  corrections  may  be 
applicable.   These  parameters  may  be  placed  anywhere  in  the  computer  and 
their  addresses  supplied  to  the  program  by  placing  them  in  the  B  registers 
in  the  following  order 

Bl  =  r  (address  of  the  crystal  radius  in  centimeters) 
B2  -  t  (address  of  the  crystal  thickness  in  centimeters) 
B3  -  h  (address  of  the  crystal  to  source  distance  in  centimeters) 
B4  =  T  (the  value  of  the  full-energy  absorption  coefficient  of  t^e 
detector  for  the  gamma  ray  in  centimeters   ) 

The  values  of  the  parameters  used  in  the  calculation  of  the  corrected 
angular  correlation  coefficients  for  the  theoretical  correlation  function 
and  the  results  of  the  calculations  are  listed  in  Table  4. 

The  program  is  available  for  future  use  on  bioctal  tape  in  machine 
language  and  on  IBM  cards  in  assembly  language.   The  bioctal  tape  contains 
all  the  necessary  subroutines  to  enable  a  complete  computations  i.e. 

1.  The  main  program  for  solution  of  Equation  (14) , 

2„  The  subroutines 
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a.  Floating  Point  Sine-Cosine 

b.  Floating  Point  Exponential 

c.  Floating  Point  Decimal  Output. 

The  results  of  the  computations  are  dumped  on  magnetic  tape  unit  (4) 
in  the  following  sequence: 


Line 

1 

xo 

Line 

2 

T2 

Line 

3 

h 

Line 

4 

I2/I0 

Line 

5 

VX0 

Line 

6 

^2 

Line 

7 

A2Q2 

Line 

8 

^4 

Line 

9 

\% 
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TABLE  4 
SOLID  ANGLE  CORRECTION  PARAMETERS 

h  (cm)  18.0 

t  (cm)  2,54 

r  (cm)  2.54 

T  (cm"1)  0.098 

I  0.0079 

I2  0.0078 

1,  0.0076 

4 

I2  0.9882 
T0 


I.  0.9609 

4 

X0 


Q2  0,9765 

Q  0.9234 

Q2A2  0.0996 

Q.A.  0.0084 

4  4 


63 


SOLID  ANGLE  CORRECTION  PROGRAM 


ORG 

12000 

REM 

BURTON  J  CONWAY  APRIL  196l 

REM 

SOLID  ANGLE  CORRECTIONS 

REM 

OF  THE  THEORETICAL  CORRELATION 

REM 

CURVE.  MUST  BE  SUPPLIED 

REM 

WITH  PARAMETER  ADDRESSES  AS 

REM 

FOLLOWS. 

REM 

Bl  ADDRESS  OF  XTAL  RADIUS. 

REM 

B2  ADDRESS  OF  XTAL  THICKMESS. 

REM 

B3  ADDRESS  OF  CRYSTAL  TO 

REM 

SOURCE  DISTANCE. 

REM 

BU  ADDRESS  OF  PARAMETER  TAU. 

HOLD 

EQU 

5500 

TEMP 

EQU 

65OO 

W0RK2 

EQU 

7500 

WORK** 

EQU 

10500 

HANS 

EQU 

11500 

SOLD 

ECU 

12500 

JAYO 

EQU 

13500 

JAY2 

EQU 

13504 

JAY*+ 

EQU 

13510 

JAY20 

ECU 

1351^ 
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L2Q00 


12001 


12002 


12003 


12004 


12006 


12007 


12010 


12011 


jayUo   equ 

13520 

EXPFL    EQU 

60600 

SINFL    ECU 

60000 

12  1  00000 

IDA 

1 

0 

20  0  12163 

STA 

0 

RAD 

12  2  00000 

LDA 

2 

0 

20  0  12164 

STA 

0 

THK 

12  3  00000 

LDA 

3 

0 

20  0  12165 

STA 

0 

HGT 

12  h   00000 

LDA 

k 

0 

20  0  12166 

ST 

0 

TAU 

12  0  12171 

LDA 

0 

ZRO 

20  5  13500 

STA 

5  JAYO 

5h   5  00030 

ISK 

5  30 

75  0  12004 

SLJ 

0 

/-l 

12  0  12171 

LDA 

0 

ZRO 

20  0  12151 

STA 

0 

DELTA 

50  1  00000 

ENI 

1 

0 

50  2  00000 

ENI 

2 

0 

50  3  00000 

ENI 

3 

0 

50  k   00000 

1 

h 

c 

50  5  00000 

ENI 

5 

0 

50  6  00000 

ENI 

6 

0 

65 


12014 


12015 


12016 


1201? 


12< 


12022 


12023 


: 


12  0  12163 
33  0  12165 
20  0  12172 
50  0  00000 
12  0  12165 
30  0  12164 
20  0  12175 
12  0  12163 
33  0  12175 
20  0  12175 
12  0  12151 
50  0  00000 

75  k  60001 
50  0  00000 
50  0  00000 
50  0  00000 
50  0  00000 
50  0  0( 
20  1  06500 
12  0 
75  k   6< 


START 


BGN 


LDA 

0 

RAD 

FDV 

0 

HGT 

STA 

0 

PLEAS 

ENI 

0 

0 

LDA 

0 

HGT 

FAD 

0  THK 

STA 

0 

ONCE 

LDA 

FDV 

0 

ONCE 

STA 

0 

ONCE 

LDA 

0 

DELTA 

ENI 

0 

0 

!  1  J 

k 

8XNFL+1 

ENI 

0 

0 

ENI 

0 

0 

1  CM 


50   0  00000 


ENI  0  0 

cjTA  i    fpRMp 

LDA 

SLJ  k  8BJFL 

ENI  0  0 

ENI  0  0 

mi 


FORMATION  OF  H/H^F 


COL 


SBi  1UE 


66 


12026 


12027 


12030 


12051 


12052 


12055 


12054 


12055 


12056 


12057 


12040 


50  0  00000 

ENI 

0 

0 

50  0  00000 

ENI 

0 

0 

20  i  05500 

STA 

1 

HOLD 

SINE  TO  STORAGE 

35  i  06500 

FDV 

1 

TEMP 

SINE/COSINE 

51  0  12175   REPL     FSB 

0 

ONCE 

SINE/COSINE -R/H4T 

22  0  12l4l 

AJP 

0 

FIRST 

JUMPS  IF  A  IS  0 

22  2  12141 

AJP 

2 

FIRST 

JUMPS  IF  A  IS  + 

51  i  ooooi 

INI 

1 

1 

12  0  12152 

IDA 

0 

DELTA+1 

30  o  12151 

FAD 

0 

DELTA 

20  0  12151 

STA 

0 

DELTA 

75  0  12017 

SU 

0 

BGN 

56  1  12105   NEXT     SIU 

1  JAY+5 

STARTS  FORMATION  1 

56  1  12075 

SIU 

1 

INCL+5 

THE  LEGENDRE  POL. 

56  1  12111 

SIU 

1  JAYT0+5 

56  1  12051 

SIU 

1 

RMK 

56  1  12074 

SIU 

1 

INCL46 

56  1  12101 

SIU 

1 

JAY-1 

56  1  12115 

SIU 

1 

JAYFR+3 

50  1  00000 

ENI 

1 

0 

12  1  06500 

LDA 

1 

TEMP 

COSB  IN  A  REG. 

52  1  06500 

FMU 

1 

TEMP 
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12041 


12042 


12043 


12044 


12045 


12046 


12047 


12050 


12051 


12052 


12053 


12054 


20  0  12173 

STA 

0 

DON 

32  0  12157 

FMU 

0 

FLT2 

31  0  12155 

FSB 

0 

FLT 

33  0  12156 

FDV 

0 

FLT1 

20  1  07500 

STA 

1 

W0RK2 

12  0  12173 

LDA 

0 

DON 

32  0  12161 

FMU 

0 

FLT4 

20  0  12174 

STA 

0 

DON+l 

12  0  12173 

LDA 

0 

DON 

32  0  12173 

FMU 

0 

DON 

32  0  12162 

FMU 

0 

FLT5 

31  0  12174 

FSB 

0 

DON+l 

30  0  12157 

FAD 

0 

FLT2 

33  0  12l6o 

FDV 

0 

FLT3 

20  1  10500 

STA 

1 

W0RK4 

50  0  00000 

ENI 

0 

0 

54  1  00000   RMK 

ISK 

1 

0 

75  0  12040 

SLJ 

0 

NEXT+4 

12  0  12164   HALF     LDA 

0 

THK 

33  1  06500 

FDV 

1 

TEMP 

20  1  11500 

STA 

1 

HANS 

50  0  00000 

ENI 

0 

0 

54  1  00000 

ISK 

1 

0 

75  0  12052 

SLJ 

0 

HALF 

P2  IN  A  REG. 


THE  A  REG. 


NUMBER  OF  ANGLES 
INSERTED  BY  PROGRAM 
FORMATION  OF  XB 


NUMBER  INSERTED 

BY  PROGRAM 
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12055  12  1  11500 

32  o  12166 

12056  20  1  11500 

50  o  ooooo 

12057  5^  l   ooooo 
75  o  12055 

12060      12  1  11500   POWER 
75  h   60600 

12061    50  0  ooooo 
50  0  ooooo 

12062  20  1  12500 

50  0  ooooo 

12063  5^  1  ooooo 
75  0  12060 

120$+      51*  2  00001 
75  0  12066 

12065  75  0  12076 

50  0  ooooo 

12066  50  1  OOOOO   INCL 

51  1  00001 

12067  12  0  12165 

33  1  06500 


LDA 

1 

HANS 

FMJ 

0 

TAU 

STA 

1 

HANS 

ENI 

0 

0 

ISK 

1 

0 

SLJ 

0 

/"2 

LDA 

1 

HANS 

SLJ 

k 

EXPFL 

ENI 

0 

0 

ENI 

0 

0 

STA 

1 

SOLD 

ENI 

0 

0 

ISK 

1 

0 

SLJ 

0 

POWER 

ISK 

2 

1 

SLJ 

0 

/+2 

SU 

0 

HANG 

ENI 

0 

0 

ENI 

1 

0 

INI 

1 

1 

LDA 

0 

HGT 

FDV 

1 

TEMP 

NUMBER  ENTERED 


STARTS  FORMATION 
OF  E  TO  THE  TAU  X 


STARTS  CALCULATION 
R  CSCB-SECB 
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12070 


12071 


12072 


12073 


1207^ 


12075 


12076 


12077 


12100 


12101 


12102 


12103 


20  1  12500 

STA 

1 

SOLD 

12  0  12163 

LDA 

0 

RAD 

33  1  05500 

FDV 

1 

HOLD 

31  1  12500 

FSB 

1 

SOLD 

32  0  12166 

FMU 

0 

TAU 

20  1  11500 

STA 

1 

HANS 

5**  i  ooooo 

ISK 

1 

0 

75  o  12067 

SLJ 

0 

/-h 

50  i  ooooo 

ENI 

1 

0 

56  1  12063 

SIU 

1 

P0WER+3 

50  1  ooooo 

ENI 

1 

0 

75  0  12060 

SU 

0 

POWER 

12  0  12155   HANG     IDA 

0 

FLT 

33  1  12500 

FDV 

1 

SOLD 

20  1  12500 

STA 

1 

SOLD 

12  0  12155 

LDA 

0 

FLT 

31  1  12500 

FSB 

1 

SOLD 

20  1  12500 

STA 

1 

SOLD 

54  1  OOOOO 

ISK 

1 

0 

75  0  12076 

SLJ 

0 

HANG 

12  0  12153   JAY 

LDA 

0 

WORK 

STARTS  CALCULATION 

32  1  12500 

FMU 

1 

SOLD 

OF  THE  INTERGAL 

32  1  05500 

FMU 

1 

HOLD 

32  0  12152 

FMU 

0 

DELTA+1 

70 


12104 


12105 


12106 


12107 


12110 


12111 


12112 


12113 


12114 


12115 


12116 


30  0  13500 

FAD 

0 

JAYO 

20  0  13500 

STA 

0 

JAYO 

54  1  00000 

ISK 

1 

0 

75  0  12102 

SLJ 

0 

/-3 

12  1  07500   JAYTO    IDA 

1 

W0RK2 

32  1  12500 

FMU 

1 

SOLD 

32  1  05500 

FMU 

1 

HOLD 

32  0  12152 

FMU 

0 

DELTA+1 

30  o  13504 

FAD 

0 

JAY2 

20  o  13504 

STA 

0 

JAY2 

54  1  00000 

ISK 

1 

0 

NUMBER  OF  TIMES 

75  0  12106 

SU 

0 

/-3 

TILL  R/H 

12  1  10500   JAYFR    LDA 

1 

W0RK4 

32  1  12500 

FMU 

1 

SOLD 

32  1  05500 

FMU 

1 

HOLD 

32  0  12152 

FMU 

0 

DELTA+1 

30  0  13510 

FAD 

0 

JAY4 

20  0  13510 

STA 

0 

JAY4 

54  1  00000 

ISK 

1 

0 

NUMBER  OF  TIMES 

75  0  12112 

SLJ 

0 

/-3 

TILL  R/H 

12  0  13504 

LDA 

0 

JAY2 

33  0  13500 

FDV 

0 

JAYO 
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12117 


12120 


12121 


1212* 


12123 


ejLX»>>  Jbt&J*' 


12125 


12126 


1212? 


12130 


12131 


12132 


20  0  1351k 

12  0  13510 
33  0  13500 
20  0  13520 
12  0  13  51k 
32  0  1351k 
20  0  1352k 
32  o  12167 

20  0  13530 
12  0  13520 

32  0  13520 
20  0  1353*+ 
32  0  12170 
20  0  135^0 
5h  1  77777 
75  0  12126 
50  0  00000 

-jh  7  420k0 
75  k  71000 
50  0  00000 

01  0  13500 

ok  0  00001 
12  0  1215k 

70  0  12131 


STA 

0  JAY20 

IDA 

0 

JAYk 

FDV 

0 

JAYO 

STA 

0 

JAYkO 

LDA 

0 

JAY20 

FMU 

0 

JAY20 

STA 

0 

JAYkO+k 

FMU 

0 

THEOR 

STA 

0 

JAYkO+10 

LDA 

0 

JAYkO 

FMU 

0 

JAYkO 

STA 

0 

JAYkO+lk 

FMU 

0 

THEOR+1 

STA 

0 

JAYkO+20 

ISK 

1 

77777 

SLJ 

0 

/ 

ENI 

0 

0 

EXF 

7 

U20k0 

SLJ 

k 

71000 

ENI 

0 

0 

01 

0 

JAYO 

Ok 

0 

1 

LDA 

0 

MIKE 

RAD 

0 

/-I 
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12133     9*  5  oooio 

75  o  12130 
1213**     54  1  77777 

75  0  12134 

12135  54  2  77777 
75  0  12135 

12136  50  o  ooooo 
74  7  42o4o 

12137  74  0  42003 
50  0  OOOOO 

i2i4o    16  o  05000 

50  o  ooooo 

12141      12  0  12147   FIRST 

20  0  12030 
121^2      12  0  12150 

20  0  12031 

12143  56  1  12054 
56  1  12057 

12144  56  1  12063 
56  1  12075 

12145  56  1  12066 

51  1  00001 


ISK 

5  10 

SLJ 

0  /-3 

ISK 

1  77777 

SLJ 

0/ 

ISK 

2  77777 

SLJ 

0/ 

ENI 

0  0 

EXF 

7  42040 

EXF 

0  42003 

ENI 

0  0 

SLS 

0  5000 

ENI 

0  0 

IDA 

0  SEC 

MODIFIES  THE  PROGRAM 

STA 

0  REPL 

WHEN  THE  VALUE  OF 

IDA 

0  SEC+1 

B     CHANGES 

STA 

0  REPL+1 

SIU 

1  HALF+2 

SIU 

1  HALF+5 

SIU 

1  POWER+3 

SIU 

1  INCL+7 

SIU 

1  INCL 

INI 

1  1 
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12146  75  0  12017 

50  0  00000 

12147  31  0  12172   SEC 
22  0  12034 

12150  22  2  12034 

51  i  ooooi 

12151  00  0  00000   DELTA 
00  0  00000 

12152  17  6  64o6i 
11  5  64570 

12153  20  0  14000   WORK 
00  0  00000 

12154  00  0  00004   MIKE 
00  0  00001 

12155  20  0  14000   FLT 
00  0  00000 

12156  20  0  24000   FLT1 
00  0  00000 

12157  20  0  26000   FLT2 
00  0  00000 

12160  20  0  44000   FLT3 
00  0  00000 

12161  20  0  57400   FLT4 
00  0  00000 


SLJ  0  BGN 

ENI  0  0 

FSB  0  PLEAS 

AJP  0  NEXT 

AJP  2  NEXT 

INI  1  1 
OCT    0 


.001 


DEC    1.0 

0  0  4 
001 

DEC  1.0 
DEC  2.0 
DEC  3.0 
DEC  8.0 
DEC    30.0 


PARAMETER  TO  MODIFY 
PROGRAM  FOR  3 
GREATER  THAN  R/H+T 
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12162 


12163 


12164 


12165 


12166 


12167 


12170 


12171 


12172 


12173 


12175 


20  0  64300 
00  0  00000 
20  0  25050 

75  3  *H217 
20  0  25050 

75  3  ^1217 
20  0  5^400 
00  0  00000 
17  7  ^6213 
20  7  12601 
17  7  46416 
25  4  02030 
17  7  14511 
51  2  62577 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 
00  0  00000 


FLT5 


RAD 


Tint 


HGT 


TAU 


THEOR 


z.yj 


PLEAS 


DON 


ONCE 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


OCT 


BSS 


BSS 


BSS 


35.0 


2.54 


2.54 


18.0 


.098 


.1020 


.00907 


0 


END 
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